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ABSTRACT 

Recently, a few peculiar Type la supernovae (SNe) that show exceptionally large peak luminosity 
have been discovered. Their luminosity requires more than 1 Mq of ^^Ni ejected during the explosion, 
suggesting that they might have originated from super-Chandrasekhar mass white dwarfs. However, 
the nature of these objects is not yet well understood. In particular, no data have been taken at 
late phases, about one year after the explosion. We report on Subaru and Keck optical spectroscopic 
and photometric observations of the SN la 2006gz, which had been classified as being one of these 
"overluminous" SNe la. The late-time behavior is distinctly different from that of normal SNe la, 
reinforcing the argument that SN 2006gz belongs to a different subclass than normal SNe la. However, 
the peculiar features found at late times are not readily connected to a large amount of ^^Ni; the SN 
is faint, and it lacks [Fc II] and [Fe III] emission. If the bulk of the radioactive energy escapes the SN 
ejecta as visual light, as is the case in normal SNe la, the mass of '^^Ni does not exceed ~ 0.3 Mq. We 
discuss several possibilities to remedy the problem. With the limited observations, however, we are 
unable to conclusively identify which process is responsible. An interesting possibility is that the bulk 
of the emission might be shifted to longer wavelengths, unlike the case in other SNe la, which might 
be related to dense C-rich regions as indicated by the early-phase data. Alternatively, it might be the 
case that SN 2006gz, though peculiar, was actually not substantially overluminous at early times. 
Subject headings: white dwarfs - radiative transfer - supernovae: individual (SN 2006gz) 



1. INTRODUCTION 

There is general agreement that Type la supernovae 
(SNe la) are thermonuclear explosions of white dwarfs 
(WDs; e.g., Nomoto, Iwamoto, & Kishimoto 1997; Hille- 
brandt & Niemeyer 2000). Thanks to the uniformity of 
their light-curve shapes after applying a correction factor 
(Phillips 1993; Phillips et al. 1999), they can be used as 
"standard candles" to measure cosmological parameters, 
which led to the discovery of the accelerating expansion 
of the universe and dark energy (Riess et al. 1998; Perl- 
mutter et al. 1999; see Filippenko 2005 for a review). 
However, the nature of their progenitor systems has not 
been resolved (e.g., Livio 2000; Hillebrandt & Niemeyer 
2000; Nomoto et al. 2003), making it difficult to rehably 
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predict potential evolutionary effects that could add sys- 
tematic errors to the determination of cosmological pa- 
rameters. 

The progenitors of normal SNe la (Branch, Fisher, & 
Nugent 1993; Li et al. 2001) are believed to be WDs 
having nearly the Chandrasekhar mass (hereafter Ch- 
SN la and Ch-WD). The recent discovery of extremely 
luminous SNe la raises the possibility that not all SNe la 
originate from a single type of progenitor system. Howell 
et al. (2006) reported that SN la 2003fg (SNLS-03D3bb) 
reached an absolute F-band magnitude of My — —19.94 
(iJo = 70 km s^i Mpc"\ f^M = 0.3, flat universe). 
Assuming that this luminosity is powered by the decay 
chain ^^Ni ^^Co ^^Fe as in other SNe la, they esti- 
mated MsBNi « 1.29 Mq (hereafter MsgNi is the mass of 
^^Ni produced and ejected during the explosion). Com- 
bining this with other elements whose existence in the 
ejecta is evident from the spectra, the ejecta and progen- 
itor masses should exceed the Chandrasekhar limit of a 
nonrotating WD. This was the first observationally based 
suggestion of an SN la from a super-Chandrasekhar WD 
(hereafter SupCh-SN la and SupCh-WD). Since then, 
two other possibly overluminous SNe la have been re- 
ported: SN 2006gz (Hicken et al. 2007) and SN 2007if 
(Yuan et al. 2007). 

All of the available data for these SNe la are only for 
the early phases, t < 100 d (hereafter t is the time since 
the explosion). At later times, SNe la enter the nebu- 
lar phase, when the Fe-rich innermost ejecta, which are 
hidden at early phases, can be directly observed (e.g., 
Axelrod 1988). In this paper, we report late-time spec- 
troscopy and photometry of the potentially overlumi- 



2 



TABLE 1 

Spectroscopy (+341 days)'' 



Filter 


Wavelengths 


Airmass 


Exposure (s) 


058 + R300 


5,800-10,200 A 


1.056 


1200 


None + B300 


4,700-9,000 A 


1.085 


1200 


Y47 + B300 


3,800-7,200 A 


1.165 


1200 



''The slit PA was deg. The flux was calibrated with the standard 
star G191B2B (Massey & Gronwall 1990). 

nous^" SN 2006gz taken with the 8.2-m Subaru telescope 
and the 10-m Keck I telescope. In §2 we present the ob- 
servations and data reduction. Results are shown in §4, 
along with some comparisons to model calculations (§3). 
SN 2006gz is intrinsically different from normal SNe la. 
However, our results are not readily interprctablc in the 
SupCh-SN la scenario, raising questions about the na- 
ture of this new subclass of objects, as discussed in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

Spectroscopy and imaging of SN 2006gz were per- 
formed on 2007 September 18 (UT dates are used 
throughout this paper) with the 8.2-m Subaru telescope 
equipped with the Faint Object Camera and Spectro- 
graph (FOCAS; Kashikawa et al. 2002). The epoch cor- 
responds to t = tpcak + 341 d, where tpcak is the time at 
S-band maximum (JD 2,454,020.2; Hicken et al. 2007) 
measured from the unknown explosion date. The field 
was also imaged on 2007 October 14 {t = tpeak + 367 
d) with the 10-m Keck I telescope equipped with the 
Low Resolution Imaging Spectrometer (LRIS; Okc et 
al. 1995). The seeing conditions were excellent on both 
nights; star profiles had full width at half-maximum in- 
tensity (FWHM) of - 0".8 and 0".7, respectively The 
same field was imaged again on 2007 November 6 by the 
Subaru/FOCAS {t = tpcak + 390 d), although the seeing 
was not good (FWHM w 1".5). 

For the Subaru spectroscopy on September 18, we took 
three spectra with exposure time 1200 s each. We used 
the O'.'8-wide slit and the R300 grism with the 058 fil- 
ter (wavelength coverage 5800 10200 A), the B300 grism 
with the Y47 filter (4700-9000 A), and the B300 grism 
with no filter (3800-7200 A), in the three separate expo- 
sures. The standard star G191B2B (Massey & Gronwall 
1990) was also observed for flux calibration. Although 
the signal-to-noise ratio (S/N) is low, we succeeded in 
obtaining a spectrum of the SN (§3). The spectroscopic 
observations are summarized in Table 1. 

Subaru/FOCAS uses an atmospheric dispersion cor- 
rector (ADC). Its performance is such that the chro- 
matic elongation due to atmospheric dispersion (Filip- 
penko 1982) is less than 0.1" within the range 3500-11000 
A at altitudes of 30-90°, so the atmospheric dispersion 
should be negligible regardless of the airmass and the 
slit angle. We confirmed that the ADC worked correctly 
during the observations; additionally, the airmass was 
low for both SN 2006gz (Table 1) and the standard star 
(~ 1.21). Thus, though the slit position angle of 0° dif- 
fered from the parallactic angle, we believe that blue vi- 
gnetting due to atmospheric dispersion was negligible. 

In this paper, we define overluminous SNe la by the peak 
luminosity, corresponding to MsgNj (mass of ^®Ni) > 1 Mq. 




Fig. 1.— A 21". 5 x 21". 5 section of the combined LRIS _R-band 
image centered on SN 2006gz taken with the 10-m Keck I telescope 
on 14 October 2007 (§2). North is up and east is to the left. The 
circle has a radius of 10 pixels, much larger than the uncertainty 

in the position of SN 2006gz (0.2 pixel). 

The total exposure time for imaging was 60s for each 
of B, V, and R on September 18, 555 s for R and 645 s for 
g on October 14, and 1000s for R on November 6. We 
obtained images of standard stars (Landolt 1992) near 
SA98-634 on September 18 and around PG0231-F051 on 
November 6 for photometric calibration. The imaging 
observations are summarized in Table 2. 

In the September Subaru images the SN was not de- 
tected. We obtained an upper limit for the SN luminosity 
as follows. First, the magnitude tuq which results in 1 
photon count per second (ADU) was derived for each 
band. Then we obtained the sky count Igky (correspond- 
ing to the sky magnitude rrisky) and the dispersion (o-gky) 
around the SN position. The Na magnitude at the de- 
tection limit is calculated as 

mum = -2.5 log ( ] + msky • (1) 

Adopting 3crsky for O'.'S x 0'.'8 binned images (i.e., the 
pixel size is nearly equal to the seeing), we estimated the 
limiting magnitude in each band as ms > 24.4, my > 
24.2, and niR > 24.0 mag. 

In October, the SN was marginally detected in both 
the g and R Keck images (Fig. 1). Since SN 2006gz 
was not readily identified in the late-time images, we 
checked the position of the SN using the stacked i?-band 
Keck image (with a total exposure time of 555 s). As- 
trometric transformation between the discovery image of 
SN 2006gz taken with the Katzman Automatic Imag- 
ing Telescope (KAIT; Filippenko et al. 2001) on 2006 
September 26 and the combined LRIS image, using 8 
stars near the location of SN 2006gz, yields a precision 
of 0.2 pixel for the position of the SN in the LRIS image. 

Figure 1 shows a 21". 5 x 21". 5 section of the LRIS 
image centered on SN 2006gz. The black circle has a 
radius of 10 pixels, which is much bigger than the as- 
trometric transformation uncertainty (0.2 pixel), and is 
used for clarity in the figure. At the center of the cir- 
cle, a faint stellar source is seen superimposed on some 
extended background emission; we regard this as a de- 
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TABLE 2 
Photometry 



UT Date (2007) Phase^ Telescope Exposure (s) Band Magnitude 



September 18 +341 Subaru 60 B > 24.4 

Subaru 60 V > 24.2 

Subaru 60 -R > 24.0 

October 14 +367 Keck 555 R 25.5 ± 0.3 

Keck 645 g 25.5 + 0.2 

November 6 +390 Subaru 1000 R > 23.3 



^Timc since S-band maximum (days). 



TABLE 3 
SN lA Models'' 



Name 




fsCE 


/56Ni 


/iME 


M(56Ni) 




SW7 


1.39 


0.18 


0.43 


0.32 


0.6 


1.40 


LW7 


1.39 


0.00 


0.72 


0.21 


1.0 


1.41 


SupCh2 


2.00 


0.21 


0.50 


0.22 


1.0 


1.60 


SupCh3 


3.00 


0.14 


0.33 


0.46 


1.0 


1.50 



^The units of mass and energy are MQand 10^^ ergs, respectively, g 
Also, /i denotes the mass fraction of different burning products (§3). 
In all of the models, the mass fraction of unburned C+O is set to be 
0.07. 



tection of SN 2006gz. Since the Subaru spectrum, al- 
though with poor S/N, shows broad features especially 
at ~7200-7300 A (§4.2), we conclude that other possi- 
bilities, such as an unrelated background object or an 
underlying H II region, are unlikely. 

The photometric reduction was performed using the 
IRAF package DAOPHOT. The SN was detected by the 
automatic star detection routine DAOFIND, and then 
the photometry was performed via PSF fitting. We ob- 
tain the R magnitude of the SN using the R magnitude 
of field stars derived from the September Subaru image. 
For the g magnitude, the Subaru B and V images for the 
field stars are used, with transformation equations from 
B and V to g given by Fukugita et al. (1996). The pho- 
tometry results in rrig = 25.5 ± 0.2 and to^ = 25.5 ± 0.3 
mag. 

In November the SN was not detected in the Subaru 
images. Because of the mediocre seeing, we did not ob- 
tain a meaningful 3a limit: mj^ > 23.3 mag. 

3. SN lA MODELS 

To quantify the observational results, we have con- 
structed four SN la models, based on the W7 model 
of Nomoto, Thielemann, & Yokoi (1984), which repro- 
duces the basic observational features of normal SNe la 
(Branch et al. 1985). Our model contains the follow- 
ing five parameters (see also Howell et al. 2006; Jeffery 
et al. 2006): M„d (progenitor WD mass), pwd (WD 
central density, set to be 3 x 10^ g cm""^ throughout 
this paper) , /ece (mass fraction of electron capture Fe- 
peak elements, e.g., ^'^Fe, ^^Fe, ^*Ni), /seNi (mass frac- 
tion of ^^Ni), and /ime (mass fraction of partially burned 
intermediate-mass elements). The kinetic energy of the 
ejecta (Ek) is given as a function of these parameters, 
as this is the nuclear energy generation reduced by the 
binding energy of the WD. We use the binding energy 
formulae from Yoon & Langer (2005), whose models in- 
clude rotating SupCh-WDs. We adopt the density dis- 
tribution of the W7 model as a function of the velocity. 
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Fig. 2. — B (blue squares), V (orange circles), and R (red tri- 
angles) light curves of SN 2006gz as compared with those of SN la 
2003du {B is the blue dashed line; V is the orange dotted line; R 
is the red solid line). Each light curve of SN 2003du is arbitrar- 
ily shifted in magnitude so that they are normalized at the peak 
and have roughly the same peak magnitude as SN 2006gz (~ 16 
mag in all bands: Hicken et al. 2007), to more easily compare 
the light-curve shapes. SN 2003du represents a typical light-curve 
shape of SNe la (e.g., SN 2001el has an almost identical light-curve 
shape; Krisciunas, et al. 2003; Stritzinger & SoUerman 2007). For 
SN 2006gz, the early-phase data (filled symbols) are from Hicken 
et al. (2007), while the late-phase date (open symbols) are from 
our Subaru (upper limit on +341 d and +390 d) and Keck (open 
symbol, +367 d) observations. SN 2003du data are from Stanishev 
et al. (2007). 

and scale it in a self-similar way: p{v) cx E-^^"^ and 

1/2 1/2 

V oc M^^ E^ . The distribution of the elements is as- 
sumed to be concentric, and ordered as follows: ECE, 
^^Ni, IME, then unburned C+O from the innermost re- 
gion. Table 3 summarizes our models: SW7 (Simplified 
W7), LW7 (Luminous W7), SupCh2 and SupCh3 (Super- 
Chandrasekhar models), characterized by different Mwd 
(progenitor WD mass), Msbni (mass of ^^Ni synthesized 
at the explosion), and Ek (kinetic energy of the expand- 
ing ejecta). 

For these models, we compute bolometric light curves 
using a one-dimensional Monte-Carlo radiation transport 
code (Cappellaro et al. 1997; Maeda et al. 2003) with 
the phenomenological opacity description for optical pho- 
tons (Kept) given by Mazzali et al. (2001a), which crudely 
takes into account the largest contribution from the Fe- 
peak elements. This description has been tested for nor- 
mal SNe la using W7-hke models (Mazzali et al. 2001a). 
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Fig. 3. — Synthetic bolometric light curve for SN la models, 
as compared with the late-time ij-band light curve of SN 2006gz 
(open triangles). Note that m/j — rriBoi = 0.4—0.5 mag for the 
models (§4.1 and Table 4). Models shown here are SW7 (thin solid; 
tpcak = 17 d), LW7 (thin dotted; tp^ak = 20 d), SupCh2 (thick 
solid; tpoak = 20 d), and SupChS (thick dotted; tp^ak = 25 d). 
Inserted figure is for the early-phase monochromatic light curves 
(Hicken et al. 2007), as compared with the same bolometric model 
light curves. We assume fi = 34.95 mag, i?(B — V^)host = 0.18 mag, 
E{B — V^)Gal = 0.063 mag, and Rv = 3.1 to convert the observed 
magnitudes to absolute magnitudes. 

The transport of 7-rays from the ^^Ni/^^Co decays is 
treated in a gray approximation with — 0.025 cm~^ 
g~^, which is very accurate for any input models (Suther- 
land & Wheeler 1984; Maeda 2006). Positron transport 
is also solved in a simplified way. We assume a phe- 
nomenological opacity Kg+ = 7 cm^ g~^, a value that ex- 
plains the slightly faster decline of late-time light curves 
of normal SNe la by increasing the fraction of escaping 
positrons (Cappellaro et al. 1997).^^ 

Synthetic nebular spectra are also computed at t = 
ipcak + 341 d (i.e., Subaru/FOCAS observation in 
September), where tpeak is obtained for each model by 
the light-curve calculations. Given the deposited lumi- 
nosity, which is obtained in the same way as in the light- 
curve calculations, ionization-recombination equilibrium 
and rate equations are solved iteratively (Mazzali et al. 
2001b; see also Maeda et al. 2006b). Since we have not 
tried to obtain detailed fits to the observed data, the 
light curve and spectrum synthesis calculations should 
be regarded only as indicative. 

4. RESULTS 

4.1. Light Curve: Rapid Fading 

Figure 2 shows the results of our photometry as com- 
bined with the early-phase light curve of Hicken et al. 
(2007). The early post-maximum decline of SN 2006gz 
was slower in all bands than that of normal SNe la 
(Hicken et al. 2007), represented here by SN 2003du. 

Sollerman et al. (2004) suggested that the light-curve be- 
havior is explained by the color evolution within the context of a 
fully trapped positron scenario. The phenomenological opacity pre- 
scription here may therefore overestimate the amount of positron 
escape. 
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Fig. 4. — Synthetic bolometric model light curves compared with 
the UBVRI bolometric light curve of SN 2006gz at early times 
(Hicken et al. 2007; see the footnote in the text). We assume 
fi = 34.95 mag, E{B - y)host = 0.18 mag, E{B - V)g!,i = 0.063 
mag, and Ry = 3.1. 

However, this is not the case at the late phase: the de- 
tection of SN 2006gz in the R band at = 25.5 ± 0.3 
mag (Keck) and the 3a upper limit in B and V (Subaru) 
show that eventually the visual luminosity of SN 2006gz 
declined more rapidly than that of other SNe la. Note 
that SN 2003du has a typical late-time light-curve shape, 
with a decline rate of 1.5-1.6 mag (100 d)~^. There have 
been several SNe la whose decline is significantly slower 
(Lair et al. 2006), opposite to the behavior seen in SN 
2006gz. 

The September Subaru photometry (+341 d), mji > 
24.0 mag, is consistent with the October Keck detection 
(+367 d) of SN 2006gz at = 25.5 ± 0.3 mag. The 
most likely magnitude at t — ipcak + 341 d is mn « 25.0, 
mv ~ 25.0, and tob « 24.8 mag^^, assuming that the 
decline rate between these two epochs follows the ^^Ni 
heating model (see Table 4 for the model prediction). 

Figure 3 shows the synthetic bolometric light curves of 
the four models as compared with the observed i?-band 
light curve. Table 4 summarizes the synthetic multi-band 
magnitudes as derived from the model spectra. The late- 
time bolometric correction is ma — msoi — 0.4-0.5 mag 
for all the models presented here. 

Figure 4 shows the synthetic bolometric light curves of 
the four models as compared with the observationally de- 
rived bolometric light curve of SN 2006gz [with Ry = 3.1, 
E{B - y)host = 0.18 mag, and E{B - V)g^x - 0.063 
mag].^^ The models with Af(^^Ni) — 1 reproduce 
the peak magnitude fairly well. The behavior before the 
peak (except SW7) is very similar for different models de- 
spite different tpcak- Thus, it seems difficult to constrain 
^pcak as accurately as tpcak = 16.6 ± 0.6 d mentioned by 
Hicken et al. (2007), and our models are consistent with 
the observed rising behavior. The different progenitor 

The B magnitude here is highly uncertain, because of the low 
S/N of the spectrum below ~4500 A. 

^•^ The observationally derived bolometric light curve is pre- 
sented at http://www.cfa.harvard.edu/supernova/SNarchive.html 
(Hicken et al. 2007). 
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TABLE 4 

Synthetic Magnitudes 



Days since B maximum Band 


SWT 


LW7 


SupCh2 


SupCh3 


Observed 


+341 


B 


24.5 


23.9 


23.7 


23.2 


> 24.4 


' 24.8)^ 




V 


23.7 


23.1 


22.8 


22.1 


> 24.2 


' 25.0)a 




R 


24.8 


24.3 


23.9 


23.0 


> 24.0 (- 


' 25.0)=- 




rriR - meol'' 


0.4 


0.4 


0.4 


0.5 






+367 


R 


25.4 


24.8 


24.4 


23.4 


25.5 






ruR - rriBol^ 


0.5 


0.4 


0.5 


0.5 







^Values in parentheses are the most probable estimates, derived as follows. For +341 d, uir is derived by extrapolating the -R-band magnitude 
at +367 d, assuming the standard ^^Ni/Co heating model. Then the spectral flux was calibrated with mR, yielding an estimate of ms and my at 
+341 d. Note that ms thus derived is highly uncertain, because of the very low S/N of the spectrum below 4,500 A. 

^Derived for the synthetic model spectra. 



masses can be distinguished after the peak. A more mas- 
sive progenitor results in a slower decline after the peak 
because of the larger diffusion time scale. From this point 
of view, Model SupCh2 with M^d = 2 Mq yields the light 
curve most consistent with the observations around the 
peak. In addition, the photospheric velocity and temper- 
ature at the peak luminosity, derived in the light-curve 
calculations, are similar in the SupCh2 and SW7 models. 
Within the uncertainties involved in the model calcula- 
tions, these are consistent with the observed character- 
istics. Details of the early-phase characteristics will be 
presented elsewhere (K. Macda, in preparation). 

At late epochs, however, the models are more than 1 
mag brighter than the observations, with the discrepancy 
reaching ^ 2 mag for SupCh2 and ~ 3 mag for SupChS 
in the V band, where the models predict strong emission 
lines (see §4.2). Indeed, only model SW7 is marginally 
consistent with the late-time photometry, while the other 
three models having Msqni = 1 Mq are too bright com- 
pared to the observations. 

The failure of the SupCh models stems from the larger 
binding energy of a WD and the higher density. The 
peak date (tpcak) measured from the explosion date and 
the peak luminosity (Lpcak) are roughly estimated by the 
following relations (e.g., Arnett 1982): 

'--^ " ( o.lem'^-0 ""Se^^^" (2) 
ipeak«7.8 X 10^3 M56Ni,©exp (^^ff^ ergs s'^ (3) 

Here; the; subscripts © and 51 mean that these values 
are expressed in units of Mq and 10^^ ergs, respectively. 
The optical depth to 7-rays from the ^^Co decay at late 
epochs {t = ttaii) is expressed as follows (e.g., Clocchiatti 
& Wheeler 1997; Maeda et al. 2003): 

T,« 1000M2,Q£-i,^(^^|i) ' . (4) 

The luminosity at ftaii (-^taii) is then estimated as 
Ltaii«1.3 X 10''^ ergs s'^Msem,© 

x(T^ + 0.035/e+)exp(^^=||^) , (5) 

where the factor 0.035 accounts for the positron contri- 
bution to the heating, with /e+ being the fraction of 
positrons trapped within the ejecta (/e+ w 1 in usual 



situations). Using these expressions and values listed in 

Table 3, the expected bolometric magnitude difference 
between the peak and the tail (ttaii = ipcak + 341 d) can 
be derived: AmBd ~ 7.1 mag for SW7 and LW7, and 
6.5 for SupCh2. This comes from the larger tpeak and 
the larger in the SupCh2 model, as a result of the 
larger mass and binding energy of a WD (i.e., larger ra- 
tios of M^^/Ek and M^^/Ek, in which Ek is reduced 
by the binding energy). This behavior, smaller AmBd in 
the SupCh models, is independent of the uncertainty in 
distance and reddening. 

Other model-independent constraints on MsgNi can be 
obtained considering the positron channel. Omitting the 
T-y term from eq. (5), we obtain an upper limit on 
-^56Ni from the requirement that this positron luminosity 
should be below the observed luminosity (Fig. 3): 

MssNi < O.3M0/-+1 . (6) 

Here we adopted mn — msoi =0.5 mag, a typical value 
in the model spectrum synthesis calculations. 

4.2. Spectrum: Missing [Fe II] and [Fe III] in the Blue 

Figure 5 shows the late-time spectrum of SN 2006gz 
and the comparison with spectra of other SNe la. In SN 
2006gz, the emissions at ~4700 A ([Fe II] AA4814, 4890, 
4905; [Fe III] AA4858, 4701) and at ~5200 A ([Fe II] 
AA5159, 5262) are extremely weak or undetected. The 
only confirmed detection is at ~7200-7300 A, which is 
interpreted as the blend of [Fe II] AA7151, 7171, 7388, 
7452 in SNe la. This feature may also be [Ca II] AA7291, 
7324 as commonly seen in core-collapse SNe (sec §5). 

Indeed, in the SUSPECT database^^ we found only 
one example of a SN la that probably shows a feature at 
~7200-7300 A as strong as the [Fe II] and [Fe III] in the 
blue. This is the underluminous SN 1991bg (Filippenko 
et al. 1992), with MggNi « 0.1 MQ(Mazzah et al. 1997). 
Thus, SN 2006gz does not appear to follow the behavior 
of normal SNe la, both in the light-curve shape (Hicken 
et al. 2007) and in the late-time spectral features. 

Figure 6 shows the synthetic spectra computed for the 
four models. There is a tendency for more-massive mod- 
els to show a weaker flux in the [Fe II]-[Fe III] emis- 
sion at <5200 A relative to the [Fe II] line near 7200 A. 
An important quantity to characterize the relative line 
strengths is the ratio of density to heating rate per mass: 
generally, as this ratio increases, ionization (competing 

http://bruford.nhn.ou.edu/ suspect /index 1. html . 
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Fig. 5. — Late-time Subaru spectrum of SN 2006gz at t = ip-|-341 
d (black line) as compared with those of SNe la 2003du (tp +377 d, 
blue line; Stanishev et al. 2007) and 1998bu (tp +329 d, green line; 
Cappellaro et al. 2001). In this comparison, the flux is normalized 
at ~7000-7500 A. The insert is the Subaru spectrum, with the flux 
(in units of 10^^** erg cm~^ A~^) renormalized to give mji = 25.0 
mag (see §4.1 and Table 4). Wavelengths are in the rest flame 
(using z = 0.02363 for SN 2006gz). The spectrum is smoothed 
with a binning width of 22.4 A. The data for SNe 2003du and 
1998bu are from the SUSPECT database (see the footnote in the 
text). 

with recombination) and temperature (competing with 
the hne emission) decrease, resulting in the stronger 
7200 A feature (corresponding to the lower ionization 
state and temperature). More-massive models have the 
larger ratio (e.g., in SupCh2 the density is a factor of 2, 
but the heating rate only a factor of 1.3, larger than in 
model SW7), but it is not sufficiently large to reproduce 
the observed, extremely large ratio of the red to the blue 
emission. 

5. DISCUSSION AND CONCLUSIONS 

We have presented late-time photometry and spec- 
troscopy of SN 2006gz obtained with the Subaru and 
Keck I telescopes. These are the first late-time observa- 
tional data for an SN la that was claimed to be overlu- 
minous at early phases. Such SNe have been suggested 
to be the explosions of SupCh-WD progenitors. 

The spectrum is characterized by the weakness of the 
[Fe II] and [Fe III] lines in the blue (< 5200 A) relative 
to emission lines in the red (either [Fe II] or [Ca II] at 
7200 A). This does not fit into the sequence of the 
late-time spectral behavior of normal SNe la, confirming 
that SN 2006gz belongs to a different subclass of objects 
than normal SNe la. Coupled with this is the problem 
of the observed faintness of SN 2006gz, especially in the 
V band. 

The SupCh2 model (with a 2 M© WD progenitor syn- 
thesizing 1 M0 of ^^Ni) is in good agreement with the 
early-phase bolometric light curve (§4.1; but see §5.1), 
although it predicts a higher late-time luminosity than 
the Ch-SN la models because of the larger M^Q-^i and 
7-ray deposition rate. The luminosity of Sup-Ch mod- 
els exceeds that observed by more than 2 mag in the V 



Fig. 6. — Synthetic late-time spectra of the SN la models. The 
model flux is converted to the observed flux using the distance 
modulus and reddening shown in the caption of Figure 3. For 
clarity, an arbitrary constant is added to the flux for each model 
(+5, +10, and +15 X 10"^^ erg cm-^ s"! A"! for LW7, SupCh2, 
and SupCh3, respectively). 

band. Furthermore, one derives M56Ni ^0.3 M©, under 
the usual assumptions that the bulk of the deposited ra- 
dioactive luminosity emerges at visual wavelengths, and 
that positrons are almost fully trapped within the ejecta. 
In what follows, we list possible solutions for this incon- 
sistency. 

5.1. Reconsidering the Super Chandrasekhar Model? 

We should first mention that the peak luminosity of 
SN 2006gz derived by Hicken et al. (2007) involves large 
uncertainty, as the host-galaxy extinction is not well con- 
strained. If the host extinction is negligible, the peak lu- 
minosity of SN 2006gz is close to that of normal SNe la 
(My = -19.19 mag). In this case, SN 2006gz might be 
a less energetic explosion of a Ch-WD, possibly because 
of insufficient burning of the C-0 layer to intermediate- 
mass elements, to be consistent with the observed slow 
light-curve evolution at early times. It is not clear if 
the early-phase spectroscopic features can be explained 
by this model: the early-phase spectra (Hicken et al. 
2007) indicate that the photospheric temperature is high 
(e.g., the weak Si II 5972) and the velocity is normal 
(z;si II « 11,000-12,000 km s~^ at B-band maximum). 
Thus, the most straightforward interpretation is that 
the peak luminosity is also high. Although the Ch-WD 
model is an interesting possibility and deserves further 
investigation, it should not affect most of our conclusions 
about the peculiar late-time behavior, as our arguments 
are based mainly on the magnitude difference between 
the peak and the tail. In what follows, we assume the 
host extinction adopted by Hicken et al. (2007), but the 
above caveat should be kept in niind.^''^ 

From the late-time data, we find that my is ~ 1.3 and 
^2.2 mag larger (that is, fainter) than expected from the 
SW7 (MgsNi = 0.6 Mq) and SupCh2 (Afgem 1-0 M©) 

For example, if the host extinction is smaller, then the esti- 
mated values for A/sgNi iu the early and late-phases both become 
smaller. 
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models, respectively. This would imply that 0.1-0.2 Mq 
of '^^Ni powers the late-time emission, smaller than the 
value expected from the light-curve peak (MsgNi ~ 
1 M0) by a factor of > 5. Such a large discrepancy can- 
not be attributed solely to the effect of viewing angle, 
even if SN 2006gz was a result of an extremely off-axis 
explosion (see, e.g., Sim ct al. 2007 for the effect of view- 
ing angle; see also Maeda, Mazzali, & Nomoto 2006a). 

Alternatively, one may hypothesize that the energy 
source at early times was not ^^Ni/Co/Fe decay. Strong 
circumstcllar interaction as seen in a peculiar class of 
SNc la (i.e., SNe la/IIa 2002ic-like events: Hamuy et al. 
2003; Deng et al. 2004; Wang et al. 2004) is not fa- 
vored, because of the lack of evidence for interaction in 
the early-phase spectra. Moreover, the SupCh models 
yield a reasonable fit to the early-time bolomet- 
ric light curve without the fine tuning of parameters. The 
early-phase spectroscopic features also seem to be con- 
sistent with SupCh models (§4.2; K. Maeda, in prepa- 
ration). Thus, ^^Ni/Co energy input is the most likely 
mechanism to power the early-phase emission. 

5.2. Positron Escape, Infrared Catastrophe, or Dust 

Formation ? 

If the early-phase emission is powered by ^ 1 M© of 
^^Ni, then one is left with the possibility that some mech- 
anism, which is not at work in normal SNe la, must be 
affecting the late-time emission and the thermal condi- 
tions within the ejecta of SN 2006gz. 

One possibility is the escape of positrons produced by 
the ^^Co decay out of the SN ejecta. The issue has 
been comprehensively explored by Milne, The, & Leis- 
ing (2001). They showed that a fraction of the positrons 
can escape for a radially combed and/or weak magnetic 
field, leading to the changing light curve at late times. 
However, at f « 400 d, this effect can account for at 
most ~ 0.5 mag, which is not large enough to completely 
remedy the present problem. 

Another possibility is the thermal catastrophe within 
the ^^Ni-rich region, which shifts the bulk of the emis- 
sion from optical to infrared wavelengths (Axelrod 1988). 
This so-called "infrared catastrophe" (IRC) is expected 
to take place after the temperature drops below a few 
1000 K, depending on the electron density. Observation- 
ally, the IRC has not been clearly detected in any SNe 
la. The IRC does not take place at i « 1 yr in any of 
our models: the heating rate per unit volume in the Fe- 
rich emission region is only a factor of ~ 1.5 smaller in 
SupCh2 than in SW7, and thus the thermal conditions 
are similar in the two models. 

Finally, it may be possible that the visual light from the 
^^Ni-rich region is converted to near-IR/mid-IR wave- 
lengths by dust formed in the C-|-0-rich region. Before 
maximum light, SN 2006gz showed the strongest C lines 
of any SN la ever observed, with little evolution in the 
absorption velocity. This indicates that a dense C+0- 
rich shell or clumps arc present in SN 2006gz, at least in 

This docs not rule out the possibility that SN 2006gz is a Ch- 
SN la with a large offset and viewed from the side of the ^^Ni blob, 
as suggested by Hillebrandt et al. (2007) for another overluminous 
SN la, SN 2003fg. Their model would have the same problem in 
reproducing the late-phase data presented in this paper, and would 
require additional processes to explain it, as do the SupCh models. 



the outermost region. The presence of carbon may be a 
key to the tmderstanding of the peculiar behavior at late 
phases, since carbon has a relatively high condensation 
temperature. 

The W7 model has n « 10^ cm~'^ at the inner edge 
of the C-l-0-rich region^® at 100 d after the explosion. 
The corresponding density in the SupCh models is a 
few times higher. Recently, Nozawa ct al. (2008) the- 
oretically investigated dust formation in SNe lb. They 
found a large amount of carbon grain formation in the 
C-enhanced outer He region, with a density comparable 
to that in the SN la models, in their model at t > 50 
d. Thus, carbon grain formation in SNe la may also be 
possible if carbon is abundant in the outermost region. 
For normal SNe la, accelerated fading like that of SN 
2006gz has never been reported. This may be consis- 
tent with the estimate that the upper limit of the carbon 
mass fraction is as low as 0.01 (e.g., Tanaka et al. 2008). 
Details of the dust formation process, addressing the dif- 
ference between SN lb and SN la models (heating rate 
by ^^Ni, different oxygen mass fraction, etc.), should be 
investigated, as well as details of expected observational 
signatures during the dust formation. In this interpreta- 
tion, a part of the outer Ca-rich region might be mixed 
with the C-hO-rich region. It may partly explain the 
relative strength of the feature at ~ 7300 A relative to 
the [Fe II] and [Fc III] in the blue, since [Ca II] arising 
from the Ca-rich region could be less diluted than the Fe 
emission, contributing to the 7300 A feature. 

5.3. Future Observations 

The possible interpretations we listed above are only 
speculative, but they predict distinct observational signa- 
tures, to be tested in future overluminous SNe la. First, 
late-time spectroscopy of other potentially overluminous 
SNe la would be extremely useful to see whether SN 
2006gz is special even among these peculiar SNe la. If 
some of them show the SN 2006gz-likc peculiarities at 
late phases, the optical light curves spanning from early 
to late times should provide a probe: (a) the IRC and 
dust scenarios predict a sudden decrease in the visual 
luminosity, (b) the positron-escape scenario would mani- 
fest itself with a relatively gentle decrease of the luminos- 
ity, and (c) other heating scenarios would not necessarily 
follow the quasi-exponential decline. 

Near-IR light curves cotild directly show the effects of 
the IRC and dust distinctly from other scenarios: we 
expect that a rapid increase in the near-IR brightness 

Khokhlov, MuUer, & Hoflich (1993) assumed that explosive 
carbon burning is ignited at the center of the sub-Chandrasekhaj: 
mass (~ 1.2 Mq) degenerate C-f-O core (with a central density as 
low as pc 2 X 10* g cm~^) surrounded by a spherical extended 
envelope. For such a configuration, they demonstrated that inter- 
action between the ejecta and the hypothesized spherical envelope 
can give rise to a dense shell. However, the WD should evolve 
until the central density becomes sufficiently high (pc fs 2 x 10^ g 
cm~^) for explosive carbon burning to occur. At this point, the 
whole mass would be in an almost hydrostatic rotating WD (Uen- 
ishi, Nomoto, & Hachisu 2003; Yoon & Langer 2004) surrounded 
by a geometrically thin accretion disk. 

This relatively dense region is formed at the interface between 
burning and non-burning layers. The density structure is almost 
frozen because the time scale of the Rayleigh- Taylor instability 
becomes longer than the time scale of the bulk expansion of the 
WD. Such a density enhancement at the (non-spherical) interface 
is thus also expected in multi-dimensional explosions. 
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should occur simultaneously with a rapid decrease in the 
visual brightness, while other scenarios do not necessarily 
predict this behavior. 

A temporal sequence of optical spectra would also be 
highly useful. In the dust formation scenario, we may 
see a transient red continuum and emission-line shifts at 
intermediate phases, as was seen in the peculiar SN lb 
2006jc (Smith, Foley, & Filippcnko 2008). 

Finally, the IRC and dust scenarios could be unam- 
biguously distinguished with near-IR through mid-IR 
spectra. Blackbody radiation from the dust particles 
should be seen in the dust scenario, while line emission 
is the dominant cooling process in the IRC scenario. 

We thank the director of Subaru, Masahiko Hayashi, 
as well as Hiroshi Terada, for kindly allowing us to ex- 
change observing time and targets in the Subaru proposal 
S07B-057 (PI: K.M.), making the SN 2006gz observations 
possible. The staffs at the Subaru and Keck Observato- 
ries are acknowledged for their excellent assistance; more- 
over, Jeff Silverman, Ryan Foley, and Maryam Modjaz 



Arnctt, W. D., 1982, ApJ, 253, 785 

Axelrod, T.S. 1988, in lAU Colloquium 108, Atmospheric 
Diagnostics in Stellar Evolution: Chemical Peculiarity, Mass 
Loss, and Explosion, (Lecture Notes in Physics, 305), ed. K. 
Nomoto (Berlin: Springer Verlag), 375 
Branch, D., Fisher, A., & Nugent, P. 1993, AJ, 106, 2383 
Branch, D., Doggett, J. B., Nomoto, K., & Thielemann, F.-K. 1985, 
ApJ, 294, 619 

Cappellaro, E., Mazzali, P. A., Benetti, S., Danziger, L J., Turatto, 

M., Delia Valle, M., & Patat, F. 1997, A&A, 328, 203 
Cappellaro, E., et al. 2001, ApJ, 549, L215 
Clocchiatti, A., & Wheeler, J. C. 1997, ApJ, 491, 375 
Deng, J., et al. 2004, ApJ, 605, L37 
FiUppenko, A. V. 1982, PASP, 94, 715 

Filippenko, A. V. 2005, in White Dwarfs: Cosmological and 
Galactic Probes, ed. E. M. Sion, S. Vennes, & H. L. Shipman 
(Dordrecht: Springer), 97 

Filippenko, A. V., Li, W., Treffers, R. R., & Modjaz, M. 2001, in 
Small- Telescope Astronomy on Global Scales, ed. W.-P. Chen, 
C. Lemme, & B. Paczynski (ASP Conf. Ser. 246; San Francisco: 
ASP), 121 

Filippenko, A. V., et al. 1992, AJ, 104, 1543 

Fukugita, M.. Icliikawa, T., Gunn, J. E., Doi, M., ShimasaJfu, K., 

& Schneider, D. P. 1996, AJ, 111, 1748 
Hamuy, M., et al. 2003, Nature, 424, 651 
Hicken, M., et al. 2007, ApJ, 669, L17 
Hillcbrandt, W., & Niemeyer, J. C. 2000, ARAA, 38, 191 
Hillcbrandt, W., Sim, S. A., & Ropkc, F. 2007, A&A, 465, 17 
Howell, D. A., et al. 2006, Nature, 443, 308 

Jeffery, D. J., Branch, D., & Baron, E. 2006, preprint (astro- 

ph/0609804) 
Kashikawa, N., et al. 2002, PASJ, 54, 819 
Khoklov, A., Miiller, E., & Hoflich, P. 1993, A&A, 270, 223 
Krisciunas, K., et al. 2003, AJ, 125, 166 

Lair, J. C, Leising, M.D., Milne, P.A., & Williams, G.G. 2006, AJ, 

132, 2024 
Landolt, A. U. 1992, AJ, 104, 340 

Li, W., Filippenko, A. V., Treffers, R. R., Riess, A. G., Hu, J., & 

Qiu, Y. 2001, ApJ, 546, 734 
Livio, M. 2000, in Type la Supernovae: Theory and Cosmology, 

ed. J. C. Niemeyer & J. W. Truran (Cambridge: Cambridge 

University Press), 33 
Maeda, K. 2006, ApJ, 644, 385 



helped with the Keck observations. We arc also grate- 
ful to Stefan Taubcnbcrgcr, Wolfgang Hillcbrandt, Fritz 
Ropke, Daniel Sauer, and Stuart Sim for useful discus- 
sions. This research is supported by World Premier In- 
ternational Research Center Initiative (WPI Initiative), 
MEXT, Japan, and by the Grant-in-Aid for Scientific 
Research of the JSPS (18104003, 18540231, 20540226, 
20840007) and MEXT (19047004, 20040004). M.T. is 
supported through the JSPS (Japan Society for the Pro- 
motion of Science) Research Fellowship for Young Scien- 
tists. A.V.F.'s supernova group at the University of Cal- 
ifornia, Berkeley is grateful for the financial support of 
National Science Foundation (NSF) grant AST-0607485 
and the TABASCO Foundation. KAIT and its ongoing 
research were made possible by donations from Sun Mi- 
crosystems, Inc., the Hewlett-Packard Company, Auto- 
Scope Corporation, Lick Observatory, the NSF, the Uni- 
versity of California, the Sylvia & Jim Katzman Founda- 
tion, and the TABASCO Foundation. This research has 
made use of the CfA Supernova Archive, which is funded 
in part by the NSF through grant AST-0606772. 



Maeda, K., Mazzali, P. A., Deng, J., Nomoto, K., Yoshii, Y., 

Tomita, H., & Kobayashi, Y. 2003, ApJ, 593, 931 
Maeda, K., MazzaU, P. A., & Nomoto, K. 2006a, ApJ, 645, 1331 
Maeda, K., Nomoto, K., Mazzali, P. A., & Deng, J. 2006b, ApJ, 
640, 854 

Massey, P., & Gronwall, C. 1990, ApJ, 358, 344 
Mazzali, P. A., et al. 1997, MNRAS, 284, 151 

Mazzali, P. A., Nomoto, K., Cappellaro, E., Nakamura, T., Umeda, 

H., & Iwamoto, K. 2001a, ApJ, 547, 988 
Mazzali, P. A., Nomoto, K., Patat, F., & Maeda, K. 2001b, ApJ, 

559, 1047 

Milne, P. A., The, L.-S., & Leising, M. D. 2001, ApJ, 559, 1019 
Nomoto, K., Thielemann, F.-K., & Yokoi, K. 1984, ApJ, 286, 644 
Nomoto, K., Iwamoto, K., & Kishimoto, N. 1997, Science, 276, 
1378 

Nomoto, K., Ucnishi, T., Kobayashi, C, Umeda, H., Ohkubo, T., 
Hachisu, L, & Kato, M. 2003, in From Twilight to Highlight: 
The Physics of Supernovae, ed. W. Hillebrandt & B. Leibundgut 
(Berlin: Springer), 115 (astro-ph/0308138) 

Nozawa, T., et al. 2008, ApJ, in press (astro-ph/0801.2015) 

Oke, J. B., et al. 1995, PASP, 107, 375 

Perlmutter, S., et al. 1999, ApJ, 517, 565 

Phillips, M. M. 1993, ApJ, 413, L105 

Phillips, M. M., Lira, P., SuntzefT, N. B., Schommer, R. A., Hamuy, 

M., & Jose, M. 1999, AJ, 118, 1766 
Riess, A. G., et al. 1998, AJ, 116, 1009 

Sim, S. A., Sauer, D. N., Ropke, F. K., & Hillebrandt, W. 2007, 

MNRAS, 378, 2 

Smith, N., Foley R. J., & Filippenko, A. V. 2008, ApJ, 680, 568 
SoUerman, J., et al. 2004, A&A, 428, 555 
Stanishcv, V., et al. 2007, A&A, 469, 645 
Stritzinger, M., & SoUerman, J. 2007, A&A, 470, LI 
Sutherland, P. G., & Wheeler, J. C. 1984, ApJ, 280, 282 
Tanaka, M., et al. 2008, ApJ, 677, 448 

Uenishi, T., Nomoto, K., & Hachisu, I. 2003, ApJ, 595, 1094 
Wang, L., Baadc, D., Hoflich, P., Wheeler, J. C, Kawabata, K., & 

Nomoto, K. 2004, ApJ, 604, L53 
Yoon, S.-C, & Langcr, N. 2004, A&A, 419, 623 
Yoon, S.-C, & Langcr, N. 2005, A&A, 435, 967 
Yuan, F., et al. 2007, ATEL 1212 



